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ORIGINAL ARTICLE 

Growth of Screen-Detected Abdominal Aortic Aneurysms 
in Men: A Bayesian Analysis 

EA Sherer 1 23 , RR Bies 13 , P Clancy 4 , PE Norman 5 and J Golledge 4 

There is considerable interindividual variability in the growth of abdominal aortic aneurysms (AAAs), but an individual's 
growth observations, risk factors, and biomarkers could potentially be used to tailor surveillance. To assess the potential for 
tailoring surveillance, this study determined the accuracy of individualized predictions of AAA size at the next surveillance 
observation. A hierarchical Bayesian model was fitted to a total of 1,732 serial ultrasound measurements from 299 men 
in whom ultrasound screening identified an AAA. The data were best described by a nonlinear model with a constant first 
derivative of the AAA growth rate with size. The area under the receiver operating characteristic (ROC) curves for predicting 
whether an AAA was >40 or >50 mm at the next observation were 0.922 and 0.979, respectively, and the median root mean 
squared error was 2.52 mm. These values were nearly identical for models with or without plasma D-dimer effects. 
CPT: Pharmacometrics & Systems Pharmacology (2012) 1, e12; doi:10.1038/psp.2012.13; advance online publication 24 October 2012 



The abdominal aorta is considered aneurysmal if its greatest 
diameter is at least 30 mm. Abdominal aortic aneurysms 
(AAAs) most commonly affect older men, with a prevalence 
rate of -2-8% in men older than 65 years. 1-3 The primary 
risk associated with AAA is aortic rupture, which is estimated 
to be fatal in 68-90% of cases. 45 AAAs are generally 
asymptomatic until rupture occurs and are mostly diagnosed 
on ultrasound or computed tomography scanning used for 
evaluation of other medical conditions or during population 
screening studies. 6-9 Because AAA rupture risk increases with 
AAA diameter, 10-12 surgical repair is generally recommended 
for AAAs measuring 50-55 mm or more depending on local 
guidelines, whereas smaller AAAs are treated conservatively 
followed by repeat imaging, usually using ultrasound. 

There is considerable inter- and intrapatient variability in the 
growth rate of AAAs with some AAAs expanding rapidly while 
others grow slowly or remain stable. 13-15 Therefore, regular 
surveillance of AAAs is recommended to monitor expansion. 
This surveillance strategy could be tailored if the growth rate of 
individual patients was more precisely known while account- 
ing for within-individual variability; surveillance could be less 
frequent for AAAs unlikely to reach 50 mm for several years 
and more frequent for those likely to reach 50 mm rapidly. One 
approach for estimating the growth rate of an individual AAA is 
to include existing observations collected during surveillance 
of that AAA to estimate its growth rate. In other words, the 
growth rate of an individual AAA is initially set to the popula- 
tion average estimate, adjusting for covariate effects, and is 
adjusted further after each of its subsequent size measure- 
ments. The objective of this work was to first develop a hier- 
archical, Bayesian model that describes growth of small AAA 
in men and then use this model to individualize AAA growth 
rate in each individual by updating the model prediction after 
each observation. The aim was to determine the accuracy of 



AAA size predictions at the next observation compared with 
the actual measurements. 

Other risk factors, biomarkers, or measures that more 
precisely predict AAA growth may be of clinical value. 16 
One example is plasma D-dimer, a factor that has been 
previously associated with AAA growth using regression 
models. 17 Although plasma D-dimer is associated with 
the population trend in AAA growth, the effect of plasma 
D-dimer (or other potential biomarkers) on individual AAA 
size predictions that are adjusted as additional surveil- 
lance size measurements is unknown. Therefore, we also 
assessed the added value of patient characteristics col- 
lected from blood tests, including plasma D-dimer, in pre- 
dicting AAA size observations. 

RESULTS 

A total of 875 men were diagnosed with a small AAA (30- 
49 mm) during the Western Australia screening study. Of 
these, 299 had both serial AAA diameter measurements and a 
D-dimer measurement and were included in this study cohort. 
In these men, the median AAA diameter at screening was 
32.7 mm (interquartile range of 30.8-36.0 mm). The men were 
followed up for a median of 5.5 years (interquartile range of 5-6 
years) and underwent a total of 1 ,732 AAA size measurements 
with a median of six (interquartile range of 6-7) measure- 
ments per patient. The median plasma D-dimer concentration 
for patients was 326 ng/nl (interquartile range of 143-786 ng/ 
ml). The demographic characteristics, risk factors at screening, 
and blood biochemistry of these men are listed in Table 1 . 

AAA growth model 

The serial AAA size measurements of individual men were 
best described by the AAA growth model with a constant first 



'Division of Clinical Pharmacology, Department of Medicine, Indiana University School of Medicine, Indianapolis, Indiana, USA; 2 Center of Excellence on Implementing 
Evidence-Based Practice, Roudebush Veterans Affairs Medical Center, Indianapolis, Indiana, USA; Indiana Clinical and Translational Science Institute, Indianapolis, 
Indiana, USA; "Vascular Biology Unit, Queensland Research Centre for Peripheral Vascular Disease, School of Medicine and Dentistry, James Cook University, 
Townsville, Australia; 5 School of Surgery, University of Western Australia, Fremantle Hospital, Fremantle, Australia. Correspondence: EA Sherer (ersherer@iupui.edu) 
Received 14 June 2012; accepted 7 September 2012; advance online publication 24 October 2012. doi:1 0.1 038/psp.201 2.1 3 



Growth of Screen-Detected Abdominal Aortic Aneurysms in Men 

Sherer eta/ 



Table 1 Characteristics of the men included in this study 
Characteristic 



Value 



Table 2 Deviance information criteria for model structures tested without 
covariate effects 



Number of men 

AAA diameter at baseline, median (q1, q3) 
Duration of follow-up, median (q1 , q3) 
AAA measurements per patient, median (q1 , q3) 
Demographic characteristics 
Age, median (q1, q3) 
Body mass index, median (q1, q3) 
Waist-to-hip ratio, median (q1, q3) 
Systolic blood pressure, median (q1, q3) 
Diastolic blood pressure, median (q1, q3) 
Medical conditions at baseline 
Diabetes, N (%) 
Hypertension, N (%) 
Coronary heart disease, N (%) 
Smoking history, N (%) 

Most cigarettes per day, median (q1 , q3) 
Pack year history, median (q1 , q3) 
Blood measurements 

D-dimer protein, median (q1 , q3) 
Glucose, median (q1, q3) 
Creatinine, median (q1, q3) 
Cholesterol, median (q1, q3) 
Triglycerides, median (q1 , q3) 
High-density lipoprotein, median (q1, q3) 
Low-density lipoprotein, median (q1, q3) 
Homocysteine, median (q1 , q3) 
C-reactive protein, median (q1, q3) 



299 

32.7 mm (30.8, 36.0) 
5.5 years (5, 6) 
6 (6, 7) 

72 years (69, 75) 
27.2 (25.1, 29.6) 
0.97 (0.93, 1.01) 
155 mm Hg (142, 169) 
91 mm Hg (83, 135) 

42 (14%) 
147 (49%) 
113(38%) 
252 (84%) 
19 (7, 30) 
29 (8, 55) 

326 ng/ml (142, 785) 
5.3 mmol/l (4.9, 5.9) 

92 umol/l (80, 111) 
4.6 mmol/l (4.0, 5.3) 
1 .2 mmol/l (0.9, 1 .7) 
1.2 mmol/l (1.1, 1.5) 

3 mmol/l (2, 3) 
13.9 umol/l (11.2, 17.7) 

2.45 mg/l (1.40, 5.06) 



AAA, abdominal aortic aneurysm. 

derivative of AAA growth rate with size (Table 2). The 
mathematical representation of this model is presented in the 
Methods section, and examples of the population and 
individual fits to data are shown in Figure 1 . In this model, the 
baseline size measurement was used to anchor the baseline 
size model parameter (10.0 mm per 10 mm increase; 95% 
credible interval (Crl), 9.95-10.01 mm). The baseline AAA 
growth rate model parameter was positively associated with 
both the baseline AAA size measurement (0.62 mm/year per 
10 mm increase; 95% Crl, 0.27-1.04 mm/year) and the 
plasma D-dimer concentration (0.36 mm/year per decade 
increase; 95% Crl, 0.04-0.65 mm/year). The first derivative of 
the AAA growth rate with size model parameter was positively 
associated with both the baseline AAA size measurement 
(0.18/year per 10 mm increase; 95% Crl, 0.03-0.31 /year per 
10 mm) and the plasma D-dimer concentration (0.15/yearper 
decade; 95% Crl, 0.05-0.25/year) but negatively associated 
with whether the patient was diabetic (-0.32/year for diabetics; 
95% Crl, -0.45 to -0.18/year). In mathematical terms, the 
expected values of the parameters in the final model for the 
baseline AAA size jgj^, baseline AAA growth rate p^ k , and 



Model 



Deviance information 
criteria (Dbar, Dhat, pD) 



Constant growth rate 7, 1 1 3 (6,573; 6,032; 540) 

Growth rate changes at a constant rate with time 6,630 (5,907; 5,185; 723) 
Growth rate changes as a linear function of time 6,509 (5,664; 4,819; 845) 
with time 

Growth rate changes as a linear function of size Did not converge 8 
with time 

Growth rate changes at a constant rate 6,487 (5,938; 5,389; 549) 

with change in size 

"The model did converge if, instead of a Wishart distribution, a series of 
inverse y distributions were used as the prior distribution for parameter 
precisions. Using these prior distributions, the model with the growth rate that 
changes at a constant rate with change in size also had the lowest deviance 
information criteria (6,532). 



r(0 



constant first derivative of AAA growth rate with size P 2J < for 
the /cth Markov chain Monte Carlo (MCMC) parameter set 
output from WinBUGS applied to the /th individual are (see 
Methods section for details on the mathematical model) 



ro.k 



Pa.k 



V,(0) 



median (V(0)) 



Y,(0) 



median(Y(0)) 
log 10 (C, (D - dimer) ) 



median(log 10 (c (D - dimer| )) 



Him A , 2,k^/ y 2,/< 



Y,(0) 



-A 



median (Y(0) 

c (D-dime,)j IOg i0 (C/ D " dimer) ) 



median(log 10 (C (D - dimer) )) 



-4 D f be,esl Diabetes,. 



where y;.(0),C/ D dimer) , and Diabetes, are the baseline AAA 
size measurement, plasma D-dimer concentration, and 
diabetes status for the /th man, respectively. The parameter 
values for the final model are given in Table 3. 

In a sensitivity analysis using categorized plasma 
D-dimer concentrations, there were significant differences 
in the baseline growth rate model parameter for patients 
with >900 ng/ml vs. those with <150 ng/ml and in the first 
derivative of the AAA growth rate with size model parameter 
for the groups with >900 and 301-900 ng/ml vs. those with 
<150 ng/ml. The deviance information criteria value for the 
model using categorized plasma D-dimer concentrations 
was nearly identical to the model using log-transformed 
continuous plasma D-dimer measurements (6,069 
categorized vs. 6,070 continuous plasma D-dimer). 

Accuracy of AAA size predictions 

A total of 1 ,677 size measurements of AAA had not been 
previously measured as >50 mm in size, and 27 men had 
an AAA that expanded to >50 mm during the study. The area 
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Figure 1 AAAs growth data (open circles), population predictions (open squares), and individual predictions for the final model (+) for 
every 10th patient by plasma D-dimer level. AAA, abdominal aortic aneurysm. 



under the receiver operating characteristic (ROC) curve for 
predicting whether AAA would be >50 mm at the next size 
measurement was 0.979 for the final model (Figure 2a). 
Similarly, a total of 1 ,228 size measurements of AAA had not 
previously been measured as >40 mm in size, and 107 men 
had an AAA that expanded to >40 mm during the study. The 
area under the ROC curve for predicting whether AAA size 
would be >40 mm at the next size measurement was 0.929 
for the final model (Figure 2b). For models without plasma 
D-dimer covariate effects, the area under the ROC curve for 
predicting whether an AAA would be >50 or >40 mm at the 
next size measurement were 0.979 and 0.922, respectively. 

The median value of the root mean square error, RMSE, of 
final model predictions for the next follow-up AAA size mea- 
surement was 2.51 mm (90% Crl of median 2.510-2.514 
mm). As shown in Table 4, the median RMSE of the final 
model predictions after a specific number of follow-up mea- 
surements ranged from 2.27 to 2.67 mm. The median RMSE 
of the model that did not include D-dimer covariate effects 
was 2.52 mm (90% Crl of median of 2.519-2.523 mm) and 
the values after a specific number of follow-up measurements 
ranged from 2.31 to 2.66 mm. 

DISCUSSION 

There is considerable inter- and intraindividual variability in 
the growth rate of small AAA, and surveillance of small AAA 



is recommended until surgical intervention is warranted. The 
surveillance interval could potentially be tailored if the growth 
rate of a particular AAA was more precisely known. We iden- 
tified an AAA growth model and found that the likelihood that 
an individual's AAA grows to >40 mm or >50 mm at the next 
surveillance measurement could be accurately predicted by 
adjusting the growth rate of an individual's AAA based on 
previous size measurements. In addition, the accuracy of the 
final model predictions for AAA size at the next surveillance 
measurement was virtually identical whether or not plasma 
D-dimer was included in the model. 

Serial AAA size measurements were included in model 
predictions to make individual predictions and because these 
measurements are available clinically as patients undergo 
AAA surveillance. If only the baseline size measurement is 
available, a classification and regression tree analysis could 
be used to aid in the identification of patients with AAA that 
are likely to grow faster or slower based on the baseline size 
measurement and plasma D-dimer concentration as previ- 
ously described. 17 In this study, it was expected that the initial 
model prediction would be identical to those of the classifi- 
cation and regression tree analysis and adjusting growth 
predictions after each subsequent measurement would 
individualize the model predictions and reduce the predic- 
tive value of the plasma D-dimer covariate effect. We found 
that the predictive value of the plasma D-dimer concentration 
was negligible. However, a prospective study during which 
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Table 3 Median values (95% credible interval) of the final model's fixed parameters, covariance matrix, and residual variability 
Fixed parameters 



/?[, y(0)) denotes covariate effect of baseline AAA size measurement on baseline AAA size model parameter 
/?J 0) denotes offset from covariate effects for baseline AAA growth rate model parameter 

^ y ' 0)) denotes covariate effect of baseline AAA size measurement on baseline AAA growth rate model parameter 

-(D-dimer) 

p\ denotes covariate effect of plasma D-dimer concentration on baseline AAA growth rate model parameter 

/?2 0) denotes offset from covariate effects for first derivative of AAA growth rate with size model parameter 

Z?^' 0 " denotes covariate effect of baseline AAA size measurement on first derivative of AAA growth rate with size model 
parameter 

-(D-dimer) 

fr 2 denotes covariate effect of plasma D-dimer concentration on first derivative of AAA growth rate with size model 

parameter 

^Diabetes) denotes covariate effect of diabetes status on first derivative of AAA growth rate with size model parameter 
Parameter covariance matrix 

O"^ = Xi 1 denotes variance of baseline AAA size model parameter 

= X 2 ,2 denotes variance of baseline AAA growth rate model parameter 

O"^ = Z 3 3 denotes variance of first derivative of AAA growth rate with size model parameter 

°« /; = 5^1,2 denotes covariance between baseline AAA size and baseline AAA growth rate model parameters 

®A>& = ^ 1 . 3 denotes covariance between baseline AAA size and first derivative of AAA growth rate with size model 
parameters 

®A.A = ^2,3 denotes covariance between baseline AAA growth rate and first derivative of AAA growth rate with size model 
parameters 

SD of the residual variability 

0.97 mm (0.93, 1.00) 



32.6 mm (32.5, 32.7) 
-1 .61 mm/year (-3.08, -0.30) 
2.03 mm/year (0.87, 3.40) 
0.90 mm/year (0.11, 1.64) 
-1 .05/year (-1 .52, -0.53) 
0.59/year (0.11, 1.03) 

0.37/year (0.13, 0.62) 

-0.32/year (-0.45, -0.18) 

0.19 mm 2 (0.13, 0.28) 
1.11 mm7year 2 (0.86, 1.42) 
0.10/year 2 (0.07, 0.13) 
0.30mm 2 /year (0.20, 0.41) 
-0.06 mm/year (-0.09, -0.03) 

-0.15 mm/year 2 (-0.22, -0.09) 



AAA, abdominal aortic aneurysm. 

D-dimer is measured at the time of repeat imaging is needed 
to better estimate the predictive effect of D-dimer in a longi- 
tudinal model. 

This model could potentially be used in the design of clini- 
cal trials for drugs that affect the growth of AAAs in men by 
simulating the likelihood of observing a significant change 
in AAA growth for a proposed trial protocol (e.g., number 
of patients needed, number and timing of observations 
per patients required, and duration of study). For example, 
we found a relatively low baseline growth rate to the mea- 
surement variability ratio (1.32 mm/year vs. 0.97 mm) and 
considerable interpatient variability in AAA growth rates. 
Qualitatively, this suggests that a clinical trial must have a 
sufficient duration and frequency of observation to separate 
the growth rate of individual subjects from the measurement 
noise and that a considerable number of patients may be 
needed to be enrolled to power the analysis. The model could 
also suggest subpopulations of patients more likely to have 
the desired AAA growth characteristics based on a patient's 
AAA size history, plasma D-dimer concentration, and diabetic 
status (assuming that the effects of the drug as a function of 



AAA size, plasma D-dimer concentration, and diabetic status 
are known). 

The dynamic AAA growth rate is consistent with the results 
from the UK Small Aneurysm Trial, which found that a qua- 
dratic model of growth rate with time best described their 
serial AAA size observation data. 15 We found that a growth 
rate that changed as the size changed, rather than with 
time, better described the Health in Men Study (HIMS) data. 
This model has the advantage of directly incorporating the 
effects of AAA size on the growth rate. Previous analysis of 
the HIMS AAA growth data using regression methods found 
no evidence of nonlinear effects. 17 This is likely because 
interindividual variability was not considered in the previous 
analysis. As shown in Figure 2, there were clearly patients 
with positive, negative, as well as no curvature. Although the 
median value of the first derivative of the AAA growth rate 
with size model parameter is small at the population level 
(0.06/year), the relatively large interindividual variability on 
this parameter (0.32/year) suggests that there are individual 
patients with significant nonlinearity in the AAA growth trajec- 
tories. This implies that the substantial differences between 
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Figure 2 Receiver operating characteristic curves for whether abdominal aortic aneurysms is observed as (a) >50 mm or (b) >40 mm for 
the first time at the next observation for models with and without D-dimer covariate effects. 



Table 4 Root mean square error of model predictions vs. data for the next 
AAA observation 



Median root mean squared error 
(90% credible interval of median value) 



observations 


Final model 


No D-dimer effects 


Overall, N= 1,732 


2.51 mm (2.510, 2.514) 


2.52 mm (2.519, 2.523) 


Individual follow-up observations 




First follow-up, 
N= 299 


2.55 mm (2.547, 2.553) 


2.57 mm (2.568, 2.574) 


Second follow-up, 
N= 285 


2.67 mm (2.64, 2.72) 


2.66 mm (2.63, 2.70) 


Third follow-up, 
W=276 


2.51 mm (2.50, 2.52) 


2.49 mm (2.48, 2.51) 


Fourth follow-up, 
W= 266 


2.51 mm (2.49, 2.54) 


2.54 mm (2.52, 2.57) 


Fifth follow-up, 
W = 250 


2.41 mm (2.39, 2.44) 


2.42 mm (2.39, 2.46) 


Sixth follow-up, 
N= 187 


2.34 mm (2.29, 2.40) 


2.33 mm (2.28, 2.41) 


Seventh follow-up, 
W=69 


2.52 mm (2.44, 2.62) 


2.51 mm (2.41, 2.63) 


Eighth follow-up, 
W = 43 


2.40 mm (2.24, 2.58) 


2.38 mm (2.22, 2.60) 



AAA, abdominal aortic aneurysm. 



individuals in this growth rate nearly averaged out on the 
population level. The data in this study are an example of 
how accounting for interpatient variability can influence 
model selection because the nonlinearity in AAA growth was 
seen at the individual level but not at the population level. We 
found that the hierarchical models with nonlinear effects had 
significantly lower deviance information criteria values than 
the linear model (7,1 13 for the linear model vs. 6,487 for the 
final model) suggesting a significantly better description of 
the data. 

Similar to previous studies, AAA growth was found to 
be positively associated with the plasma D-dimer concen- 
tration 17 and AAA diameter at screening 1415 but negatively 
associated with diabetes. 17-20 Similar to the meta-analysis 
of Sweeting et a/., 18 AAA growth was independent of age, 
sex, and mean arterial pressure. In contrast to the findings 



of Sweeting et a/., 18 we found no effect of smoking on 
AAA growth. 

A log-transformation of plasma D-dimer concentration had 
the best association with growth rate model parameters sug- 
gesting that differences in AAA growth were predominantly 
driven by patients with the highest plasma D-dimer concentra- 
tions. This finding was reenforced by a sensitivity analysis in 
which plasma D-dimer values were categorized. This analysis 
demonstrated differences in the baseline growth rate model 
parameter only between patients with the lowest (<150 ng/ml) 
and highest (>900 ng/ml) plasma D-dimer concentrations. The 
baseline AAA size measurement was also positively associ- 
ated with the change in growth rate as well as the baseline 
growth rate. This result is internally consistent; because the 
baseline growth rate parameter was higher for larger AAAs, it 
follows that the growth rate would increase as size increased. 
Finally, diabetes was only associated with the first derivative 
of the growth rate with AAA size model parameter but not the 
baseline growth rate model parameter. In a previous study, 
Vega de Ceniga et al. 20 found that diabetes was associated 
with slower AAA expansion only for AAAs with an initial size 
of 4-4.9 cm. It is possible that the effect of diabetes becomes 
more evident as AAAs increase in size which was difficult to 
be assessed in this study because most of the AAAs were <4 
cm. In patients with diabetes, Golledge et al." found a stan- 
dardized regression coefficient of -0.17 between the overall 
growth rate and the diabetic status of the patient. We found no 
associations between growth rate and smoking or hyperten- 
sion that have previously been reported as significant. 1518 

The population baseline growth rate model parameter of 
1.32 mm/year was slower than the 2.60 mm/year and 2.81 
mm/year growth rates reported in the UK Small Aneurysm 
Trial 15 and the Chichester trial, 14 respectively. This difference 
is most likely due to the smaller initial AAA sizes in the HIMS 
study relative to those in the other studies (32.7 vs. 43 mm of 
median baseline AAA size in HIMS vs. UK Small Aneurysm 
Trial, respectively). Adjusting for initial AAA size, the predicted 
initial growth rate of an AAA that was 43 mm at baseline is 
1 .98 mm/year, but this remains lower than the initial growth 
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rate reported in the UK Small Aneurysm Trial and Chichester 
studies. The reasons for these disparities remain unknown 
but could relate to differences in the population or methods 
of identifying subjects. Men included in this study were identi- 
fied through screening rather than recruitment and they had 
to return for another blood sampling years after their initial 
diagnosis; it is possible that the cohort examined in this study 
is focused on more slowly growing AAAs, because those that 
were initially larger at diagnosis or grew faster would have 
already required AAA repair. 

This study has several limitations. Foremost is that the 
D-dimer measurement was not collected at baseline and 
only measured at one time point. Therefore, the precise 
value of D-dimer in predicting AAA growth is not clear from 
this study; we found that D-dimer was indicative of AAA that 
did grow and not necessarily an AAA that will grow. In addi- 
tion, while AAA growth rate was associated with D-dimer 
level, it is possible that the plasma D-dimer was associ- 
ated with the larger size rather than the growth rate. A post 
hoc analysis found a small, positive correlation between the 
final AAA size observations and the D-dimer measurements 
(Pearson's r= 0.26). Longitudinal measurements of D-dimer 
are needed to determine whether D-dimer varies with AAA 
size or instantaneous growth rate. It is possible that longi- 
tudinal measures of plasma D-dimer could be predictive of 
instantaneous AAA growth. The study was also conducted 
in a single geographic location with a predominantly white 
population, on males only, and on small, screen-detected 
AAAs. The generalizability of the findings to other patient 
populations and to non-screen detected or large AAA 
is uncertain. The predictive ability of the model was also 
tested on surveillance intervals of mostly 0.5- or 1-year 
intervals and the accuracy for longer intervals is unknown. 
Finally, because of the relatively small sample size of the 
cohort, the accuracy of the predictions was tested on the 
same sample used for model development. The accuracy of 
the model when applied to an external cohort is unknown. 

METHODS 

Study setting and participants. The cohort of men included 
in this study consists of individuals who had a small (30-49 
mm) AAA identified during the HIMS screening for AAA; had 
a follow-up ultrasound scan for surveillance of AAA; and pro- 
vided a blood sample during the HIMS follow-up survey which 
was previously used to measure plasma D-dimer. Details of 
the design of the HIMS screening study and follow-up survey 
have previously been reported. 172122 Briefly, men included in 
the HIMS were originally part of a population screening study 
for AAA in Perth, Western Australia. The study was con- 
ducted from 1996 to 1999 and involved men of age 65-83 
years. These men received an ultrasound for AAA screen- 
ing and completed a history and lifestyle survey. Men were 
followed up as part of the HIMS between 2001 and 2004. 
Surviving men from the screening cohort were invited to pro- 
vide a blood sample. The Human Research Ethics Commit- 
tee at the University of Western Australia approved the ethics 
of the HIMS, and all subjects provided written and informed 
consent before participating in the study. The data analysis 



protocol for this study was approved by the Indiana University 
Institutional Review Board. 

AAA imaging. AAA size was measured as the greatest 
diameter of the infrarenal aorta using ultrasound (Toshiba 
Capasee with a 3.75 MHz probe, Toshiba Australia, North 
Ryde, Australia). The reproducibility of ultrasound measure- 
ments and interobserver variability were assessed every 
4 months by re-imaging randomly selected patients. As 
previously reported, there were no significant differences 
between observers, and 95% of the measurement differ- 
ences were <3 mm. 7 AAA surveillance via ultrasound imag- 
ing was recommended every 6 months for patients with 
AAA > 40 mm in diameter or every 12 months for patients 
with AAA diameter of 30-39 mm. Referral for consideration 
of surgical intervention was recommended for patients with 
AAA > 50 mm. 

Risk factors. At the time of AAA identification, all men com- 
pleted a questionnaire designed to collect relevant patient 
history and lifestyle factors. The questionnaire assessed self- 
reported smoking status and number of cigarettes; history of, 
or treatment for, coronary heart disease, peripheral arterial 
disease, hypertension, stroke, dyslipidemia, and diabetes; 
and family history of AAA (first-degree relatives only). Age, 
height, weight, blood pressure, and the circumference of hips 
and waist of subjects were also measured. 

Blood assays. Blood was collected between 2001 and 2004, 
and serum and plasma were isolated and stored at -80°C 
Plasma D-dimer was measured by enzyme-linked immu- 
nosorbent assay as previously described. 17 The interassay 
coefficient of variation for plasma D-dimer concentrations 
was 2-3%. Serum glucose, creatinine, total cholesterol, trig- 
lycerides, high-density lipoprotein, low-density lipoprotein, 
homocysteine, and C-reactive protein were measured by 
automated assays as previously reported. 2324 

Model of AAA growth. A hierarchical Bayesian model 
with individual-level growth parameters was fitted to the 
serial AAA size data using WinBUGS 1.4.3 25 software 
(see Supplementary Data online) called from R using 
the R2WinBUGS package. 26 Five chains were used for 
each analysis with each chain having a burn-in period of 
10,000 iterations and postburn-in period of 20,000, keeping 
every 10th sample for a total sample of 10,000 iterations. 
Gelman-Rubin diagnostics were used to confirm conver- 
gence of the model (see Supplementary Figures S1 and 
S2 online). 

The base (or covariate free) model was selected by 
comparing the deviance information criteria and parameter 
Crl among several potential model structures. The candidates 
included models with a constant growth rate with time; a 
growth rate whose first derivative with time was a constant 
rate, a linear function of time, and a linear function of AAA 
size; and a growth rate whose first derivative with size was 
a constant rate. The following differential equations give the 
expected AAA diameter, Y,,k(t), for the klh MCMC iteration 
applied to the /th individual as a function of time. The predicted 
AAA diameter at this time, Yi,k(f), is assumed to be normally 
distributed about the expected value where e~N{o,dfj is 
the amount of deviation. 



CPT: Pharmacometrics & Systems Pharmacology 



Growth of Screen-Detected Abdominal Aortic Aneurysms in Men 

Sherer eta/ 



Constant AAA growth rate. Individual-level baseline size, 
fl%\, and growth rate, parameters 



^^'l, V ; (0) = ^» 



df 



or 



First derivative of the AAA growth rate with time is constant. 
Individual-level baseline size, fl^ k , baseline growth rate, fly,, 
and constant first derivative of AAA growth rate with time, flf , 
parameters 

^ = *;,(*), Y,AQ) = fi% 

] -^ k , r,,(0) = ^> 



or 



dr 



Y,At) = flZ+fl^t+/^f + e 



First derivative of AAA growth rate with time is a linear 
function of time. Individual-level baseline size, flly, baseline 
growth rate, fly\; and constant, flf^, and linear, fl^ k , rates of 
change in AAA growth rate with time 



^ = ^(0, n,(o)=^> 

df 



2/^ + 6/?<<>f, U0) = fl\' k 



or 



First derivative of the AAA growth rate with time is a linear 
function of AAA size. Individual-level baseline size, fl^ k ; 
baseline growth rate, fl\\ \ and constant, flf M , and a positive 
linear with size, j3f k , rates of change in AAA growth rate with 
time 



6Yi, k (t) 
df 

d^(0 _ 
df 



Y'A0) = flZ 



'■PZ+fl%Ym, r tJl (0) = fll 



or 



Y»(t): 



JflZ~+ 



L+1 



exp 



exp 



First derivative of AAA growth rate with size is constant. 
Individual-level baseline size, flf k , baseline AAA growth 
rate, flf^ k , and constant first derivative of AAA growth rate 
with size, fll] k . 



dx 



l J 2,k' 



r u (/u(0)) =/?<;> 



or 



0,k 



9(0 
\k 



?(0 A 

0,/t 



^k- 



f J 2,k"o,k 



'2,k 



The regression parameters for an individual patient 
P*' = \flo ) k'fl\k'---'fl l nk\ are selected from a multivariate 
normal distribution p^' ~ Ay^p^xj where p^' is the vector of 
expected parameter values for /cth MCMC parameter set to the 

covariate effects of the /'th individual and 2 is the covariance 
matrix. Noninformative normal distributions were used as 
prior distributions for the parameters, a Wishart distribution 
for parameter precisions, and an inverse gamma distribution 
for the precision of the normally distributed residual error. 

After selection of the base model structure, covariate 
effects on model parameters were tested using a stepwise 
forward addition and backward elimination method. Additive, 
proportional, and power functions were tested for each 
continuous covariate and additive and proportional functions 
were tested for each discrete covariate (see Table 1 for patient 
demographic, medical conditions, and blood measurement 
covariates tested). In addition, plasma D-dimer was also 
tested as a discrete covariate by grouping according to 
concentration (<150, 151-300, 301-900, and >900 ng/ml). 
Covariate effects whose 95% Crls did not cross 0 in single 
covariate analyses were included in an intermediate model. 
Covariate effects whose 95% Crls did not cross 0 in the 
intermediate model were eliminated and additional iterations 
of the forward addition and backward elimination were 
performed until there were no changes. 

When the baseline size measurement was included as a 
covariate on model parameters, this measurement was not 
used in the Bayesian analysis. In other words, the baseline 
size measurement was used to adjust model parameters either 
through covariate effects or in a Bayesian manner but not both. 

Accuracy of AAA size predictions 

To quantify the predictive capability of the final model, we 
calculated the RMSE and the area under the ROC curve of 
the model-predicted likelihood that an AAA was at least 40 mm 
(after which the surveillance interval was reduced from 12 to 
6 months) or 50 mm (after which referral for consideration of 
surgical intervention was recommended) at the next observation. 
To assess the added value of D-dimer for predicting AAA size 
observations, we compared these outcomes for the final model, 
which includes patient-specific D-dimer measurements, vs. 
those for a model without plasma D-dimer effects. 
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The RMSE compared AAA size measurements — where Y.. 

and t. were the observed AAA size truncated to the nearest 

u 

0.1 mm and amount of time after the screening observation, 
respectively, for the y'th observation of the /th individual, A/ 
was the number of AAA size observations of the /th individual, 
and x was the vector of covariate measurements for the /th 
individual — with the corresponding model-predicted size 
truncated to the nearest 0.1 mm, YiiT tne MCMC 

parameter set, and P t , output from WinBUGS. The model 
predictions for y'th observation were conditional on all previous 
size observations because the likelihood of a particular 
parameter set was conditional on the prior observations 



RMSE = 



L2, 9 I^X!°=rPr(Ain/-i >V.. ;,x. 

K 



-(0.1 mm) 

V,m \t;M..x 



v(0.1mm)i 



A particular value of Yij,k |f,,;P^;x ( was randomly 
sampled based on the approximation of MCMC iterations of 
the posterior predictive distribution 



v(0.1 mm) i 



: Pr(y < Yij, k \t u ; p t ,x,. < y + 0.1 mm 



where yis in 0.1 mm increments beginning at 0 and 



Pr 



\y\y y. 



n^Pr 



-10.1mm) 

y*Y v * =y v \ 

t,j] $ k ;x, < y + 0.1 mm 



io,ooo££f°njrj,Pr 



—{0.1 mm) 

ysv v * =Y,,\ 
f y ; p„; x, < y + 0.1 mm 



The RMSE calculation was replicated 1 ,000 times to deter- 
mine a 90% Crl about the median value. The RMSE for a spe- 
cific observation number (i.e., value of y) was also retained if 
there were at least 40 observations. 

The ROCcurvewasgenerated based on clinicalobservations 
for two outcomes used during the HIMS: the likelihood of a 
first size observation of at least 40 mm and the likelihood of 
a first size observation of at least 50 mm. After a particular 
value of Yij was selected, this value was compared with both 
40 and 50 mm. For the 40 mm threshold, the model-predicted 
probability that the AAA size at a particular observation was 
>40 mm was the fraction of simulations for that observation 
that was >40 mm. To generate the ROC curve, the threshold 
above which the model-predicted probability of being >40 
mm was defined as "predicting" this size was varied and the 
actual observation of whether AAA was at least 40 mm at 
the next observation with the model prediction. Observations 
for patients with AAA that had already crossed 40 mm were 
excluded from that analysis. An analogous approach was 
used for the 50 mm threshold. 
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Study Highlights 

WHAT ISTHE CURRENT KNOWLEDGE ON THE 
TOPIC? 

There is considerable interindividual variability in AAA 
growth, and plasma D-dimer concentration is associated 
with growth. 

WHAT QUESTION DIDTHIS STUDY ADDRESS? 

This study sought to develop a hierarchical, Bayesian 
model that describes growth of small AAA and then to 
determine whether the model can accurately predict AAA 
size at the next observation. 

WHATTHIS STUDY ADDSTO OUR KNOWLEDGE 

A nonlinear growth model for AAA in which the growth 
rate changes as AAA grow or shrink best fitted the data. 
This model has advantages of being consistent with the 
baseline AAA size covariate effect and not relying on 
temporal growth rate changes. The model accurately 
predicted whether an AAA would be >40 or >50 mm at the 
next surveillance observation. 

HOWTHIS MIGHT CHANGE CLINICAL 
PHARMACOLOGY ANDTHERAPEUTICS 

This model could potentially be used to tailor surveillance 
intervals based on model-predicted risk. We found that, 
although plasma D-dimer was associated with AAA 
growth, there was little added value of a single plasma 
D-dimer concentration for predicting AAA size. 
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